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Abstract. The question of the K factor in transversely polarised Drell-Yan (DY) processes is examined. 
The transverse-spin case is peculiar for the absence of a reference point in deeply inelastic scattering (DIS). 
Therefore, in order to study more fully the possible effects of higher-order corrections on DY asymmetries, 
a DIS definition for transversity is devised using a hypothetical scalar (Higgs-like) vertex. The results show 
that some care may be required in interpreting experimentally extracted partonic transversity, in particular 
when comparing with model calculations or predictions. 



1 Introduction 

The theoretical framework for describing transversity (at 
the basic level of partonic processes, QCD evolution, ra- 
diative effects etc.) is now solid [1] and a number of ex- 
periments aimed at its measurement are on-line or under 
development: HERMES [2], COMPASS [3] and the RHIC 
spin programme [4]; there are also proposals for Drell- 
Yan (DY) measurements with polarised antiprotons in the 
High Energy Storage Ring at GSI [5, 6] (related prelimi- 
nary theoretical studies have been made regarding access 
to transversity in J/tp production [7,8]). 

Transversity is the last remaining piece in the partonic 
jig-saw puzzle composing the hadronic picture. However, 
the standard procedure of adopting deeply inelastic scat- 
tering (DIS) as the process to define parton densities at 
the next-to-leading order (NLO) cannot be extended to 
transversity in a simple manner since it does not con- 
tribute to DIS. Furthermore, transverse-spin effects are 
notoriously surprising; e.g., see the large and (historically) 
unexpected single-spin asymmetries [9,10]. Such consid- 
erations render imperative the complete understanding 
of NLO corrections in DY before attempts are made to 
extract the partonic transversity distributions. See, e.g., 
Ref. [11] for a detailed discussion of transversity and also 
single-spin asymmetries. 

One might instead consider double-spin asymmetries 
Att for other processes, such as: p^p^ — > jet+X, j+X etc. 
Unfortunately, however, predictions for Att always turn 
out to be very small [12, 13], so that measuring transver- 
sity directly appears feasible only in doubly polarised pp 
interactions. 

Since all QCD and electroweak vertices conserve quark 
chirality, transversity actually decouples from DIS. Chiral- 



ity flip is not a problem though if the quark lines connect 
to different hadrons as in, e.g., the DY process. Unfortu- 
nately, there is a caveat to accessing transversity in DY: 
Hikasa's theorem [14], which states that, owing to chi- 
ral symmetry, transversity effects vanish upon integrating 
over the lepton-pair azimuth. No simple proof of the the- 
orem exists (it has to do with the 7-matrix algebra). Let 
us now make a few observations based on these properties 
of transversity: 

1. the only "gold-plated" process in which transversity 
may be measured directly (i.e., without the need of 
more-or-less exotic fragmentation functions) is DY; 

2. Hikasa's theorem implies the use of a slightly less than 
fully inclusive process, in as much as one angle must 
be left unintegrated; 

3. in the case of transversity asymmetries, helicity conser- 
vation may not necessarily provide the usual safeguard 
against large K factors. 

The above have non-trivial implications with respect to 
the measurement of transversity in DY and interpretation 
of the results. 



1. In the absence of a DIS reference point, there is no im- 
mediate way of fully evaluating the possible importance of 
higher-order QCD corrections. The K factors are known 
to be large at the level of cross-sections in both the un- 
polarised [15] and helicity-dependent [16] cases. However, 
in the helicity case the large corrections cancel in the 
asymmetry [16]. To a large extent this cancellation can 
be traced to the conservation of helicity along fermion 
lines in gauge theories — the 0(cv s ) Wilson coefficient for 
the DY process is identical for the helicity-dependent and 
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-independent pieces, as too are the leading order (LO) 
anomalous dimensions. Note that in the case of heavy- 
flavour production the corrections to the helicity asym- 
metry are large, precisely because mass terms introduce 
helicity flip, destroying the usual protection. 



2 + 3. The coefficient function for transversely polarised 
DY differs significantly from the other two cases [17, 18]. 
Moreover, the LO anomalous dimensions differ — there is 
no corresponding conserved quantity or sum-rule. 

Given the marked differences from the other two cases, 
it may be useful to examine the question of DY K factors 
for transversity. In order to do this, it is clearly necessary 
to find some suitable DIS-like process as a reference point. 
The principal requirement is a spin-flip mechanism. There 
are two obvious possibilities a priori: either a quark mass 
term or a scalar vertex. Now, of course, DIS with trans- 
versely polarised leptons and nucleons should be consid- 
ered and therefore the twist-three structure function 172 
for general reviews) is the natural object of study (see 
e.g., Rcfs. [19,20]. It turns out, however, that although 
transversity is intimately related to the evolution of 52 
(the relevant operator is indeed proportional to the quark 
mass [21-23]), at the level of direct contribution to po- 
larised DIS it actually cancels against other higher-twist 
contributions owing to the equations of motion, see for 
example Rcf. [24] . Although the calculation is rather del- 
icate, the possibility of defining a coefficient function for 
transversity via its role in the evolution of gi has been ex- 
amined [25], with similar results to those presented here. 

A simpler and more direct approach is to identify a 
DIS-like process in which a scalar particle plays a role. 
Since the Higgs boson does indeed interact with quarks 
(as with leptons too), the obvious solution to the prob- 
lem is a gedanken process in which the exchange is no 
longer via the electroweak gauge fields but via the Higgs 
particle. To be precise, in order to obtain the required 
single spin-flip, Higgs- Vector interference diagrams actu- 
ally need to be considered. Of course, there is no intended 
suggestion here that such a process should really be mea- 
sured, but merely that it forms a suitable basis for a theo- 
retical cross-check. We should remark that such a process 
has effectively already been exploited for the calculation 
of h\(x) itself [26], on the basis of a suggestion by Jaffe. 
In any case, various tests will be performed to ensure that 
the results do not depend on the specific nature of the 
vertex introduced. 

Before moving on to the calculation, let us spend a few 
more words on the physical significance of the isT-factor. 
While at a theoretical level the meaning of higher-order 
corrections to any given process is clear, at a phenomeno- 
logical level in the parton-model there is an inherent am- 
biguity owing to the necessary input of the parton densi- 
ties themselves. Indeed, the if-factor was used historically 
to represent the discrepancy between experimental results 
for the DY cross-section and the LO theoretical predic- 
tions based on parton densities extracted from DIS and 
thus by definition (as in Rcfs. [15, 16]) the phenomeno- 
logical i-T-factor is the translation factor from DIS to DY. 



Therefore, since all model calculations or estimates of par- 
tonic transversity densities rely to some extent on DIS for 
overall normalisation or determination of model param- 
eters, self-consistency would require a procedure of the 
type to be described here. Since, furthermore, the overall 
if-factor so-defined receives large contributions from both 
DIS and DY, this is a non-trivial observation. The pecu- 
liar structure of transversity leaves room for very different 
corrections as compared to the spin-averaged or helicity- 
dependent cases, for both DIS and DY independently. 

Thus, in the following section the calculations are de- 
scribed, the Higgs-Vector interference mechanism is ex- 
amined in detail and NLO calculation of the related Wil- 
son coefficients is performed. The known results for the 
DY process are discussed and finally the relevant K fac- 
tors are extracted. In the closing section some conclusions 
are drawn and comments relevant to future measurements 
of transversity via DY scattering are made. 1 



2 The Calculation 

2.1 Drell-Yan cross-section and asymmetries 

It is now standard to define the helicity- and transversity- 
weighted cross-sections by 
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where the prefixes A and At indicate longitudinal-spin 
(or helicity) and transverse-spin (or transversity) depen- 
dence respectively, ± refer to initial-state proton helici- 
ties and |,l to transverse polarisations. The double-spin 
asymmetries are then 



and 
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The large NLO corrections afflict both the numerators and 
denominators. The question is to what extent they are 
correlated, i.e., to what extent they are the same and 
thus cancel in the ratio. 

Turning then to the calculation of the K factor, the 
procedure will be essentially identical to that followed in 
earlier work [15, 16] and thus we shall not dwell on the 
general technicalities, save for those points that are sig- 
nificantly different in the case of transverse polarisation. 
The first peculiar aspect to be exploited is that, owing to 



1 Owing to correction of an error in the code used for numer- 
ical estimates, the results shown here are a little less dramatic 
than those presented by the author in past conferences. 
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the charge-conjugation properties of the relevant operator, 
the evolution of transversity is of the flavour non-singlet 
(NS) type. In the NS case the effect of higher-order correc- 
tions may be represented in the following schematic way: 



F(x,t) 



a s (Q 
2vr 



a s (Q 2 ) 
2tt 



9/(y 5 *)> ( 3 ) 



where t = ln(Q 2 / fi 2 ), with Q 2 the virtuality of the photon, 
Qf(y,t) and Qf are respectively the parton density and 
charge of quark flavour /, P is the universal quark-quark 
splitting function and C the process-dependent Wilson co- 
efficient. The quantity F(x,t) on the left-hand side then 
represents a generic (flavour NS) structure function and 
the three terms inside the square brackets on the right- 
hand side represent: 1. the LO point-like contribution; 
2. the leading-logarithmic correction; and 3. the NLO cor- 
rection. It is this last that is of interest here. 

To NLO the DY cross-section for pp scattering is ex- 
pressed in term of parton densities as follows: 
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where r = Q 2 / 's, Q 2 is the invariant mass squared of the 
lepton pair and s is the total hadron centre-of-mass en- 
ergy squared. In Eq. @| 1^2 are the momentum fractions 
carried by the (anti)quarks inside hadrons 1 and 2 respec- 
tively. It is then the difference between the DIS correc- 
tions, with which the NLO parton distributions are de- 
fined via Eq. (J2J), and NLO DY corrections that consti- 
tutes the phenomenological K factor. 

The leading-logarithmic splitting functions P are well 
known [27-32] and may be expressed in the following com- 
pact form: 



AP(z) = P(z) = C* F 



and 



1 + z 2 



1-z 



A T P(z)=P(z)-C F (l-z), 



(5a) 
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The definition of the so-called "plus" regularisation is re- 
called in Appendix [X] Already then it is evident that 
although fermion-helicity conservation guarantees identi- 
cal evolution for NS spin-averaged and helicity-weightcd 
quark densities, the same does not hold for the transver- 
sity case. 

The problem now is to calculate the coefficient C(z) 
of the third term in Eq. @. This must be done for both 
the DIS and DY processes. As is well-known, a large part 
of the DY K factor can be attributed to the change from 




Fig. 1. The DIS "handbag" diagram for a photon-Higgs in- 
terference process. 



a space-like Q 2 in DIS to time-like in DY. However, this 
is not the only origin of large corrections and one should 
be concerned that the transversity case, with the extra 
requirement on the final-state phase space in DY, might 
introduce important differences. 



2.2 The Drell-Yan Process 

Since the results for the DY process are known [17, 18], 
it is perhaps better to begin with this coefficient. The 
virtual photon decays into a final lepton pair, of which 
then the azimuthal angle must be left unintegratcd in the 
case of transversity. In the modified minimal subtraction 
(MS) scheme the results for the unpolarised [15], helicity 
[16] and transversity [18] DY coefficient functions are as 
follows: 2 
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It is important to note that although the large 7r 2 terms 
and indeed the coefficient of the (5-function (which indeed 
constitute the bulk of the K factor) appear invariant here, 



6z In" 1 



in the transversity case there is a new term, — (i- z ) 
found in the others. 



not 



2.3 Deeply Inelastic Scattering 

In order to accommodate spin-flip in the standard DIS 
"handbag" diagram, one of the vertices should involve a 
Higgs-like scalar, see Fig. ^ The contribution of this di- 



Here and in what follows only the unpolarised result will be 
presented in full while the helicity and transversity cases will 
be shown as differences with respect to the unpolarised case. 
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Fig. 2. The two diagram types contributing to the NLO DIS 
hard partonic q-y*(H*) — > gq scattering subprocess, the dotted 
line represents either a virtual photon or Higgs. 



agram can be expressed in terms of a structure that will 
be called hi here for brevity. Projecting with j5p$ T then 
leads to 



W» = h 1 (x,Q 2 ) 



p-q 



(7) 



2.3.1 Real-gluon contributions 



The NLO Wilson coefficient may now be calculated from 
the diagrams in Fig. [3 The process to be calculated is, of 
course, still photon-Higgs interference. The use of dimen- 
sional regularisation poses the problem of dealing with 
75, which naturally arises in the case of polarised DIS 
(for both hclicity and transversity) owing to the projec- 
tor j5p^ T . Note that for the qq — > j*g DY subprocess 
this is not a problem since both the quark and antiquark 
bring one power of 75, which then cancels before calcu- 
lating any traces. The technique adopted here is that of 
defining a fully anti-commuting but non-cyclic 75, see for 
example Ref. [33]; see also Ref. [34] for a detailed discus- 
sion of this technique. Consistency then requires that all 
traces be evaluated from the same reference point, which 
in the usual DIS case is unambiguously one of the photon 
vertices. Here the scalar vertex could also be chosen — an 
explicit check shows that there is no ambiguity in the re- 
sults. 

In the MS scheme (see Appendix for a working defi- 
nition), adopting the above-mentioned 75 scheme and sup- 
pressing (for clarity) a common factor 



C, 



r(l - 2e) 



on the right-hand side of all equations, the results for the 
real contributions of the diagrams in Fig. [21 are: 
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where the C ms ~ R (z) are defined to be the combined quan- 
tities 

C ms (z) = tP{z) + C ms {z) 1 (9) 

with P(z) and C(z) being replaced respectively by AP(z) 
and Cl(z) etc., where necessary. In the first equation, for 
C(z), an additional contribution 3z due to Fl has been 
included to give the correction corresponding to the use of 
F2 to define q(x) [15]. Moreover, in Eq. (IHcfl the remain- 
ing e is due to the difference in splitting functions and 
disappears in the final expression for the full coefficient. 

To extract the desired coefficients C DIS_R (z), the vir- 
tual corrections must now, of course, also be added. First 
however, note that the results for the unpolarised and he- 
licity cases agree with previous calculations, [15] and [16] 
respectively. Note also that the results for the various 
cases are (not surprisingly) similar: while the coefficient 
for hi is a little different (owing to the finite residues of 
the ultraviolet divergences, which lead to different split- 
ting functions), the infrared (IR) double poles in e are 
identical (and in any case cancel with the virtual contri- 
butions) and the single poles themselves are, of course, to 
be absorbed into the scale-dependent parton densities. In 
particular, there is no trace of the 6 * ln „? term found in 
the DY coefficient. 



(l-z) 



2.3.2 Virtual-gluon contributions 

The virtual contributions can be partially gleaned from 
the literature; however, the scalar-vertex correction re- 
mains to be evaluated and this requires a some care. The 
real and virtual contributions are separately gauge invari- 
ant and a natural choice (as in Refs. [15, 16] and other 
cited work) is the Landau gauge, where it is only the ver- 
tex correction that need be calculated. Schematically, 



and 
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where 1 represents the "bare" scalar vertex. In MS the 
following results are then obtained: 



S v = C F 
and 
S S = C F 
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Noting that these corrections multiply 6(1 — z), one imme- 
diately sees that the double poles 1/e 2 , of IR origin, cancel 
against the real diagrams, just as they should. However, 
note also that the single-pole structure is manifestly dif- 
ferent; let us examine this a little more closely. 

There is a substantial difference between a vector and a 
scalar vertex: the former is related to a conserved current, 
the latter not. Thus, the vector current is not renormalised 
while the scalar does receive radiative corrections. In other 
words, one must also take into account the renormalisation 
of the coupling constant (i.e., the quark mass in the true 
Higgs case) associated with the vertices in consideration. 3 
Indeed, the simplest way to evaluate the contribution is 
to calculate the quark- mass renormalisation, including the 
constant pieces. In the MS scheme the standard calcula- 
tion gives 



a s (Q 2 ) 
2tt 



1 



r 2 (l - e) 3 

^(1-26) e l-2e' 



(12) 



Including this as a contribution to the virtual corrections, 
one finally obtains 

4 uU = S V . (13) 

Thus, combining real and virtual contributions, the com- 
plete set of coefficients for DIS are 
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2.4 The K-Factor Results 

The DY and DIS coefficients can now be combined to pro- 
vide a theoretical K factor. Note that the DIS coefficient 
appears with a factor two in the required difference; this 



3 This observation was made by Bliimlein [35] in regard of 
similar calculations aimed at evaluating the evolution kernel. 
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Fig. 3. The DY-DIS difference in the Mellin moments of the 
Wilson coefficients (i.e., the K factor) for the three leading- 
twist densities. 



merely reflects the two quarks (or rather quark-antiquark) 
in the initial state for DY. The final results are 4 



C UY (z)-2C ms (z) = C F 2(l + z 2 ) 
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where the origins of the large differences in the last line 
may thus be traced in part to the different phase-space 
restrictions in the transversity case and in part to the 
residues due to the different splitting functions. It is per- 
haps worth reminding the reader that the bulk of the 
large K factor in the unpolarised and helicity cases (the 
7r 2 terms) comes from the necessary continuation of Q 2 
from space-like (in DIS) to time- like (in DY). However, 
the transversity correction contains other non-negligible 
pieces. 

For a first visual comparison, Fig. |21 shows the Mellin 
moments, defined by 



(«) 



dx x n f(x), 



(16) 



of the above differences in the Wilson coefficients between 
DY and DIS for the three leading-twist densities. While it 
is clear that the difference is generally only a little larger 
for transversity, the rapidly growing difference as n — ► 
(equivalent in z space to z — > 0) with respect to both 
the unpolarised and helicity cases is particularly striking. 



4 A brief review of the results presented here may be found 
in a recent contributed talk [36]. 
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Note also that the small- z (small-n) region counts heavily 
in the convolution integrals, therefore the corrections to 
the transversity- weighted cross-section can be significantly 
larger. 

To provide an idea of the effect these corrections might 
have on an experimentally measured DY asymmetry, in 
Fig-EJthe helicity and transversity asymmetries for purely 
NS contributions (in both the numerator and denomi- 
nator) are plotted as functions of r. The asymmetries 
are shown for a centre-of-mass energy yfs — 200 GeV, 
corresponding to recent polarised RHIC energies [4]. For 
the transversity distribution we have taken a conservative 
starting point of Arq(x, Qq) = Aq(x,Qq). The evolution 
of the distributions has then been performed to LO here 
as the coefficient differences are scheme independent and 
the effect of higher orders on the K factor is negligible. 
The size of the shift due to the K factor is two to three 
times as large in the case of transversity with respect to 
helicity and typically reaches values of the order of 10% 
for these energies. It should be noted, however, that there 
is an automatic limitation of the if -factor difference {e.g., 
with respect to large values of a s for low energies) owing 
to the presence of the it 2 terms; when the if -factor differ- 
ence between numerator and denominator becomes large 
in absolute terms so too do the overall K factors them- 
selves, which to some extent cancels or dilutes the effect. 
Very similar corrections to those presented are obtained 
for energies corresponding to the proposed experiments at 
GSI. Repeating the calculations for various centre-of-mass 
energies, one finds that the corrections are maximal at 
y/s « 40 GeV, where for the transversity case they reach 
as much as 15%. 

2.5 Cross-checks 

A couple of simple cross-checks may be made on the in- 
fluence of the scalar vertex itself. First of all, there is now 
the possibility of a purely Higgs, unpolarised DIS process 
(i.e., in which both vertices are scalar). The contribution 
of the real diagrams is 

C Dls - K > s {z) = C DIS " R (z) - 2 - 3z, (17) 

which, combined with the virtual corrections already dis- 
cussed, gives 

C» IS > s ( z ) = C ms (z)-2-3z. (18) 

The new coefficient in Eq. (|18f) could be used in place 
of the usual unpolarised correction to define the parton 
distributions. 

Secondly, there is also similarly a possible purely scalar 
DY-like process; the NLO correction to the unpolarised 
cross-section in this case is found to be 

C DY ' S (z) = C DY (z) + C F 2(l-z). (19) 

Moreover, Hikasa's theorem is avoided here owing to the 
presence of the scalar vertices and a transverse-spin asym- 
metry is present even after integrating over the lepton- 
pair azimuthal angle. The NLO correction to the scalar 



transversity asymmetry is 

C^ Y ' S (z) = C DY (z) 

+ C F (l-z)[4-41n(l-z) + 21nz]. (20) 

In none of the above cases does the scalar vertex intro- 
duce large correction differences with respect to the vec- 
tor. Indeed, in the last case of a purely scalar DY process 
(both spin averaged and transversely polarised), the only 
differences are residues of the difference in the LO split- 
ting functions, as indicated by the form and the overall 
factor (1 — z). 



3 Conclusions 

In order to appraise the real nature of the DY K factor 
in the case of transversity asymmetries, we have examined 
gedanken processes involving scalar vertices. This allows a 
natural DIS definition for the partonic densities Axq(x). 
Such a definition allows a direct connection with model 
estimates based on knowledge of parton densities derived 
essential from precisely DIS. Typical examples might be 
models in which at some low Q 2 scale transversity- and 
hclicity-weighted densities are naturally equal or others in 
which the Soffer bound [37] is found to be saturated, again 
at some low scale. In all such cases the spin-averaged and 
hclicity-weighted densities used to set the starting point 
are obviously and naturally taken directly from DIS. 

The results presented here provide a measure of the 
reliability of model predictions, without, of course, repre- 
senting a rigorous estimate, in as much as the reference 
processes are partially fictitious and in any case are not 
precisely those normally adopted. However, we have seen 
that in general the corrections are not excessively large 
although they may be significantly larger than in the he- 
licity case. Moreover, comparison of the Mellin moments 
indicates that in kinematical configurations in which low 
z = r jx\Xi dominates there could be very important cor- 
rections. On the other hand, many of the differences be- 
tween NLO coefficients vanish numerically for z — ► 1 and 
so safe kinematical configurations certainly exist. 

In closing then, although apparently fairly well un- 
der control, the question of NLO perturbative corrections 
in the case of transversely polarised DY processes clearly 
deserves more study, in particular, where the kinematics 
might be such experimentally as to favour the dangerous 
low-z region. 



Acknowledgments 

Loop calculations were performed with the aid of FORM 
version 3.1 [38]. The author would like to thank Johannes 
Blumlein for very useful conversations in the past concern- 
ing the calculation of transversity anomalous dimensions 
via photon-scalar interference diagrams. 



Philip G. Ratcliffe: Transversity K Factors for Drell-Yan Processes 



7 



0.4 



0.3 



0.2 



0.1 













LO 






NLO 


0.0001 


0.001 


j 0.01 



(a) 



0.4 



0.3 



0.2 



0.1 









LO 




NLO 



0.0001 



0.001 

(b) 



0.01 



Fig. 4. The doubly polarised pp helicity and transversity asymmetries All and Att for purely valence-driven DY at LO and 
NLO as functions of r = Q 2 /s, for = 200 GeV. 



A Plus-regularised distributions 

The so-called "plus-regularised" distributions are defined 
via integrals with a smooth test function /(y): 



dyf(y) 



g(y) 



= / dy 



i-y 



g{y) 



/(l)fdA (21) 
Jo 1-V 



where g(y) is well-behaved as y — > 1. 

B Modified minimal subtraction scheme: 
implementation 

The minimal subtraction (MS) scheme is defined, in con- 
junction with dimensional regularisation, as the removal 
of all simple poles in l/e (double and higher poles due to 
IR divergences are cancelled automatically between real 
and virtual contributions). However, common residual fi- 
nite contributions are always left due to the appearance 
of the factor (47r) e .T(e). Expanding to O(e), one obtains 



(47r) e r(e) ~ r(l + e) [1 + e 1h(4tt)] 

1 

e 



+ ln(47r) — 7e- 



(22) 



The MS scheme then augments MS by subtracting the two 
e-independent terms above. 

Thus, MS may be implemented to 0(a s ) by defining 
the Feynman virtual momentum- integral measure [d n fc] to 
include a factor l/r(l + e) and by removing a factor (47r) e . 
In other words, the plain MS definition may be substituted 
with the following: 



[d n k] 



d n k 



(47r) e r(l + e) J (2tt) 



(23) 



Consequently, the definition of the phase-space integral for 
a final two-body state (as in jq — > qg) must be modified 
analogously to 



1 

8^ 



dy [y(l - y)] 



(24) 



where y = ^(l + cosfl) and 9 is the partonic centre-of-mass 
scattering angle (in this case between the incoming q and 
outgoing g). 
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